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Abstract 

We provide an exhaustive classification of self-dual four-dimensional grav- 
itational instantons foliated with three-dimensional homogeneous spaces, i.e. 
homogeneous self-dual metrics on four-dimensional Euclidean spaces admit- 
ting a Bianchi simply transitive isometry group. The classification pattern 
is based on the algebra homomorphisms relating the Bianchi group and the 
duality group 50(3). New and general solutions are found for Bianchi III. 
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1 Introduction 



Classical solutions of general relativity in diverse dimensions have been analyzed 
over the years from various perspectives. Instanton-like configurations have in par- 
ticular attracted attention because of their potential role in the determination of 
transition amplitudes in quantum gravity. Similarly, their real-time counterpart 
turns out to be useful in the Hartle-Hawking formulation of quantum cosmology, 
even though general relativity would ultimately require a ultra-violet completion, 
possibly provided by strings. In fact, several classical solutions turn out to be em- 
beddable in string theory, sometimes even as exact backgrounds, and instantons 
have certainly played a major role in addressing many issues, among which are 
those related to (super) symmetry breaking. 

Following the paradigm of self-dual Yang-Mills instantons [TJ, self-duality has 
been successfully implemented in four-dimensional general relativity. In order to 
be operational, self-duality must be accompanied by some specific ansatz for the 
geometry M.. The usual ansatz is to assume M. topologically RxE 3 and, further, 
the leaves S 3 to be homogeneous spaces, admitting at least a three-dimensional 
group of motions G3. 

The first solution obtained according to the above pattern is due to Eguchi 
and Hanson [21 [3] under the assumption of Bianchi IX geometry i.e. with SU(2) 
isometry. This gravitational instanton was an alternative to the earlier Taub-NUT 
metric, constructed from a different perspective though jl]. Both were self-dual 
and S*?7(2)-homogeneous, with isometry enhancement to SU(2) x U(l). General 
extensions of Eguchi-Hanson solution to strict SU(2) were obtained soon after [5], 
observing that the equations were a special case of the Euler top - with some pe- 
culiar inertia momenta - known as the Lagrange system [6]. A similar achievement 
for the Taub-NUT instanton revealed far more involved, and a particular solution 
was finally obtained in Ref. [7]. The difficulty to solve in full generality the cor- 
responding equations was later understood in terms of non-algebraic integrability 
properties, as it was realized [8] that the system at hand had already been set 
by Darboux [9] , and solved extensively by Halphen [TO} ITT] , more than a century 
before, in terms of modular formal 

The above account for the Bianchi IX group raises immediately two questions: 
do other Bianchi groups possess similar solutions and what is the classification 
principle behind the appearance of distinguished classes of equations such as La- 

*A translation of the general Halphcn-Darboux solutions in terms of gravitational instantons 
can be found in [12] , They are all plagued with naked singularities, except, marginally, for the 
particular solution of Atiyah-Hitchin [7] . 
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grange versus Darboux-Halphen systems? Despite the large amount of information 
accumulated so far (see e.g. [13] HI]) and the physical interest of some - even sim- 
ple - solutions like Kasner's Bianchi I [15], a precise and definite answer to these 
questions was still missing. The aim of the present work is to tame the plethora of 
scattered results under a simple classification principle. Our analysis is performed 
along the lines announced in [TH]. It is general and exhaustive, and makes no 
assumptions on the geometry, other than those already quoted above. In partic- 
ular, the issue of the choice of a diagonal versus non-diagonal metric in a given 
G3-invariant frame is treated with care, as opposed to some former, more cavalier 
approaches for non-unimodular Bianchi groups. As a bonus, this enables us to 
discover a new solution for Bianchi III, completing thereby the existing landscape. 

Let us summarize the method and the results. Self-dual vacuum solutions sat- 
isfy f2 = tl where f2 is the Riemann curvature two-fornjl. These are second-order 
equations and are equivalent to the first-order set obtained with the connection 
one-form: u = Cj + A. The one-form A appears as a "constant of motion" of the 
second-order system. It stands as the anti-self-dual so (3) part of the Levi-Civita 
connection and must be flat. The program is thus cast as follows: (i) find all 
possible flat so (3) connections over G3, and (ii) for each of them, write the corre- 
sponding first-order equations, and find the most general solution. The latter can 
represent a bona fide geometry, but in most cases it is spoiled by naked singular- 
ities, if not everywhere degenerate like in all non-unimodular groups, except for 
Bianchi III. 

In Sec. [2] we provide some technical tools, useful to set our philosophy for the 
subsequent developments. Section [3] contains the core of the classification principle 
(point (i) above), whereas the exhaustive solution search (point (ii)) is presented 
in Sec. HI A last section ([5]) is devoted to the subtle issue of rotating (locally) the 
frame where u = Cj + A into a frame where the connection is genuinely self-dual. 
Conclusions follow and a summary of all G3-invariant metrics finally available is 
presented in the appendix. 



2 Self-duality can be imposed alternatively on the Weyl tensor, and leads thus to solutions 
with cosmological constant like Fubini-Study or Pedersen [T7] . Note also that anti-self-solutions 
are obtained by parity or time reversal. 
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2 Notations and general considerations 

Inspired by applications to homogeneous cosmology, we consider Euclidean spaces 
admitting a G3 group of motion acting transitiveljH on 3-dimensional invariant 
subspaces of the form M = 1 x S 3 . The metric on these spaces can always be 
locally written as (see for instance Ref. [18], page 516, for a proof): 

d s 2 = N(t) 2 dt 2 + s ab e a e b . (1) 

Let us introduce the co-frame defined by 0° = N(t) dt and a = 0^(t) a a , where 
<j a are invariant 1— forms: 

1 



do a = ^y No\ (2) 



From the latter we obtain 



d<d a = 6 a a dtAa a + l -6 a a c a M a p A a 7 (3) 

= AQ b + \e a a c%e^2Q b aq c (4) 

providing the non vanishing structure coefficients: 

p a ob = ^« = -/n , (5) 

0% c = e a a c a ^e^2 = -p a cb , (6) 

and connection coefficients 

lOab = -(PabO + Pbao) = lOba, (7) 

labO = ^(PabO — Pbao) = — 7&a0> (8) 

labc = -{fiabc + fibca ~ ficab) = —Jbac- (9) 



With respect to the basis {dt, a a } the metric components are g 00 = N(t) 2 , g 0a = 

ianb 



and g a s = 5 a b9®9 b . The Levi-Civita connection one- forms are: 



QUij — TUijo dt -\- TUij a (J C 



3 The three-dimensional group G3 acts simply transitively on the leaves, endowed thus with 
the structure of a group manifold. Hence we exclude H-^^H^ x S 1 or S 2 x S 1 , which are the 
alternatives to the nine Bianchi classes. 
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Finally, we need the expression of the curvature components, which are: 



dw 0i + ujQk A u\ 

Wm n dt Aa a + \w 0ia c a /?7 a p A a 1 



— U7 0ka tJ k i0 dt A a a + G7ofc/3G7 fc j 7 <7^ A <T 7 , 

Vt jk = dojjk + oJjo A uP k + ojji A u l k 



(10) 



= rijjfca rft A a a + -Wjfca c ^ a 13 A a 7 + rujoa^^cr A <x^ 

H-iUj^o^'fea dt A a a — zuji a w l k0 dt A cr a + Wji a zu l k/3 a a A a 13 . (11) 

With respect to this basis, the anti-self dual components of the curvature 2- 
form read as: 



jk ( • l l \ 

ttOia — WQka^iQ ~ n e i \ W jka + Z&jioW ka — TUji Q ZU fcQ J 



1 

— ( 

2 



dt A o c 



+- 



ol k k j k 

ttOi a C + WQ k pW j 7 — WQ k ^W — 



-zu 



jkaC^ + WjopW kl + WjipW kl 



(12) 



a 13 A a 7 



or, introducing 



_ jk 



(13) 



(14) 



Let us emphasize however that the one-forms 9 k are not completely fixed, but de- 
fined up to time-dependent 0(3) transformations and time-independent GL(3, R) 
transformations, we have indeed not yet fixed the basis of invariant forms. 



3 Self-duality equations 

Self-duality requires that 

fl<H - \^ji0 k = 0. (15) 

We see that to trivially obtain first integrals from this equation a sufficient (and 
necessary) condition is to require that 

Wjfeo = ia 9 k — 9 ka 9f — 9i a — —g (X ^9 i , (16) 
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after which (TH|) becomes equivalent to 

4 = I ia c a p 7 + 4 k I jP I kl = (17) 

The gauge freedom allows us to always satisfy this condition. Indeed if we have 
a solution 9 aa , using a time dependent rotation we obtain 6 aa := O a b (t)9b a which 
will fulfill the symmetry condition if 

6a=\o:{6 a c e b a -e ca e ab ). (is) 

This is a kinematic problem aiming to determine a sequence of rotations knowing 
the angular velocity so that it always admits a solution. Then the only trans- 
formations still possible correspond to the choice of the initial conditions of the 
kinematical problem, denoted from now on by the matrix O h al and of arbitrary 
GL(3, R) time independent transformations, denoted from now on by the matrix 

A/3, 

0aa H- O a %^i ■ (19) 

In the gauge defined by (1161) . the first integrals (Ii a = const.) furnish the first order 
equations we need to solve: 

:= -I ia (20) 

where 6 = \J det(g and the symmetric tensor density n 7/i = and vector 

cip are defined by the relations 

c a M e^ = 2{n aiM + e a ^ap). (21) 

In the following we shall call the t coordinate time and speak about evolution, to 
describe this flow, though the framework is Euclidean. 



N 



Z $ t fx $ ^ , )'~\ 



72 



7/x _ 



a n e 



Q\ia, 



-S 1 
2 



Using time re-parameterizations, we may choose a specific N without loss of 
generality. Equations ff20l strongly suggest to adopt, in a first step at least, the 
gauge N = 9. Let us notice that in order to preserve this gauge condition when 
we make a transformation (fTU|) we have to simultaneously rescale t by a factor 
1/ det(A^); then the time variable may still only be changed by a shift 1 1— >■ t + to- 
But before trying to integrate equations (TSUI) it is mandatory to solve the constraint 
equations: 

I ia c a M + ef 7^4 7 = . (22) 

The solutions of these equations describe homomorphisms from the Lie algebra of 
the homogeneity group g 3 into the Lie algebra 50 (3): 

I : 03 ->so(3). (23) 
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Depending on the subalgebra of so (3) on which we project g3 different simplifica- 
tions can occur. Let us remind that the only subalgebras of so (3) are the trivial 
ones: so (3) and the null one {0} and the one-dimensional ones : R (all equivalent, 
i.e. linked by internal conjugation). Thus if the constants of motion I ia are not 
identically zero (in which case we have the trivial homomorphism mapping g 3 onto 
{0}) either g 3 = so(3) or the algebra g3 has a 2 dimensional ideal. 

These remarks lead to the following Bianchi types to be considered, according 
to the rank of the matrix I ia : 

• rank 3 (maximal) type IX , 

• rank 2 impossible, 

• rank 1 types I, II, III, IV, V,VI,VII, 

• rank all Bianchi types . 

But it remains to examine if all these cases can effectively be obtained i.e. if there 
are no obstructions coming from (Tl6|) . Let us notice that using (|20|) . condition 
( TTO|) can be written as 

hJ* ~ IjJi + 2 a p 9 pk e kij = . (24) 

From this equation we see that for rank 0, we must have a p = 0, which is the 
statement that rank solutions must be from Bianchi class A, i.e. types I, II, 
VIo, VIIo, VIII and IX, see table CD for our conventions. For the rank 1 cases, 
as shown in the next section, the only solutions whose metric determinant does 
not everywhere vanish are those of Bianchi class A and type Illfl Thus if we also 
require the metric determinant to not everywhere vanish, the above list is reduced 
to 

• rank 3 (maximal) type IX , 

• rank 2 impossible, 

• rank 1 types I, II, III, VI ,VII , 

• rank types I, II, VI , VII , VIII, XL 

Thus, according to the rank of the mapping I we obtain the classification of solu- 
tions in the next section. 

4 To see this, first decompose Z; a as Ii a = X Ja where I a and \ are two triples. For all types of 
Bianchi class B, except type III, the only solution to the constraint is to take I a proportional 
to a a . In this case, one may show that (|24[) requires the determinant of Of to vanish. 
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4 All self-dual solutions 

In this section, we obtain all the self-dual solutions categorized by the rank of /. 
Some of them have been discussed in the past in [TBI 114^1 . On a similar ground 
of simplicity, the only rank 3 solution coming from the type IX Bianchi algebra 
has also been studied in [TBI fT4"] and various other places in the past. We finally 
deal with the new cases, solutions of rank 1. Thanks to residual symmetries 
and constraints coming from the self-duality requirement, we simplify as most 
as possible the equations and present the most general solutions as well as some 
particular simplified solutions when available. 

Rank zero All are zero. The Bianchi type must be of class A, (a p = 0), 
i.e. types I, II, VI , VII , VIII and IX. Equations (12TJI) imply that the symmetry 
condition is always satisfied : 



At an initial time, one may simultaneously diagonalize g^ v and put in its canon- 
ical form (see tabled]). Thus at initial time the co- frame components 6*j 7 may be 
chosen diagonal, and they remain so during their evolution. The explicit integra- 
tion of the evolution equations has been described a long time ago in Ref. [T3| fT4"]. 
here this class of self-dual spaces is presented in appendix A for completeness. 

Rank three Of course only Bianchi IX is possible for this case of maximal rank. 
In an appropriate frame the matrix n is the identity matrix. The self-duality condi- 
tion f[2"2l can be reinterpreted as a relation between three 3-dimensional Euclidean 
vectors : 



and two similar relations obtained by circular permutation of the indices 1,2,3. 
Thus these three vectors are perpendicular to each other and actually define an 
orthonormal frame of E 3 . Acting with an appropriate element of 0(3), on the flat 
indices i, we may assume that — {ha) is the identity matrix. The integration of 
(1201) . under these assumptions, was first discussed in [6j [9j [TTJ [TT] and the solution 
is given in appendix A. 

5 However in Ref. [T31 Q3] the metrics were assumed to be a priori diagonal, a technical 
assumption just introduced to facilitate the integration of the flow equations. Here we justify 
that they furnish the most general ones in case of rank cases, and when necessary shall discuss 
the integration of the relevant non-diagonal metrics. 




(25) 



h = I 2 x I 3 



(26) 



S 



Table 1: Canonical structure constants for the different Bianchi groups 



Type 


a 


n 1 


n 2 


n 3 


Usual name 


Class A 


I 














Translations 


II 





I 








Galilean 


VIIo 





I 


I 





Euclidean 


VI 





1 


-1 





Poincare 


IX 





1 


1 


1 


Rotation 


VIII 





1 


1 


— I 


Lorentz 


Class B 


V 


1 













IV 


1 


1 










VII, 


h > 


1 


1 







VI^! 


h > 


1 


-1 







Ill EE VIi 


I 


I 


— I 








Rank one In that case there exists two triplets I a and Aj such that 7j a = Aj/ Q . 
Condition (|22p reduces to 

I«<% = 0. (27) 

Bianchi class A In this case, condition (|27[) is equivalent to 

I a n a P = 0, (28) 

so the matrix (n al3 ) has to be singular. Consequently the Bianchi type has to be 
I, II, VIo or VIIo- The symmetry condition (|24|) becomes then: 

\ iPj -\ jPi = with P k = 6%I a . (29) 

It is easy to see that if the initial conditions satisfy this condition, then the equation 
will automatically be preserved throughout the evolution. This can be seen by 
computing the evolution of pi using (1201) to obtain: 

Pi = {\n^g^ u + \ k p k ) Pi . (30) 
Thus we find no restrictions on the existence of solutions to our self-dual problem. 
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In order to find the rank 1 solutions for Bianchi class A, let us start with the 
most general expression of the matrix of frame components : 

(31) 

It can be simplified by making use of the transformations (1T9]) . Using an 0(3) 
transformation we may assume A = (0,0,1) and, without altering the canonical 
value of the structure constants, using GL(3,M.) transformations set 

• I a = (0, 0, 1) for Bianchi 1,11 

• I a = (0, 0, Z) for Bianchi VI ,VII 

Moreover, for Bianchi type I, II, VIo and VIIo (Bianchi type of class A), without 
loss of generality, we may also assume the matrix (l3"Tj) diagonal. This can be done 
using the transformations (1191) . which now leave I a and Aj invariant as well as 
the canonical values of the structure constants. The integration of the resulting 
self-duality equations was given in ref [T3| [H] . 

For illustrative purpose let us discuss the Bianchi type VIo- The integration 
proceeds as follows. First we observe that ( E5J) imposes u — v — 0. Equations f[3"0"j) 
then insure that if u, v are chosen zero at the initial time, they will remain null. 
Then we remark that the evolution equations imply that the functions q and s are 
proportional to c: 

q{t) = Qc(t), s{t) = Sc(t) (32) 

and that the transformations ( TT9"j) of the frame that preserves the canonical choice 
of the structure constants, of I a and A are given by: 

/ ±cos(0) sin(0) 0\ / a 7 \ 

0= Tsin(0) cos(0) and A 1 = I (3 a 5 \ . (33) 

\ 1 / \ 001 / 

Transformations ( 133]) allow to simplify the frame components, without loss of 
generality. First we may put Q = S = i.e. q(t) = s(t) = by performing a 
transformation with O = Id and a = 1, /3 = 0, 7 = — Q and 5 = —S. Then we 
still have the freedom to make transformations (f3"3"|) with 7 = 5 = 0. Performing, 
at an arbitrary time to, a new transformation with still O = Id and a = p cos (<£>), 
j3 = p sin ((f) where p > is arbitrary and (p determined bjjf) 

• ,„ \ 2(q r +po b ) 

8 " 1(2 ^ = -WTWTJTJ) m 



6 We use the subscript to indicate the value of the corresponding function at the instant 
considered : a (to) — a o 5 etc. 
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we obtain a new frame such that now aor + b p = 0. Let us notice that (134)) 
makes sense because from (do ±^o) 2 + (fro ±Po) 2 > we deduce that its right-hand 
side member is always between —1 and +1. Then choosing A = Id and O with 
such that 

sin(0) p + sin(0) b = = — sin(0) a + cos(0) r (35) 
we have obtained a diagonal frame at time to- 

It remains then five functions to be determined: a, b, r, p and c. These functions 
have to satisfy the equations: 

aA-2Kb 

a = — > (36) 
■ ftA + 2/ta 

& = ^-5 > ( 37 ) 



2k 2 
pA + 2nr 

r A — 2np 
2^2 ' 



P = —2 > ^ 



(39) 



A-2Zk . . 



with 



A = [b 2 + r 2 ) - (a 2 +p 2 ), and K = ab-pr. (41) 

It is obvious from these equations that k is a constant. Moreover, the fact that we 
may diagonalize the frame at t = to implies that it remains diagonal during all of 
its evolution: thus r(t) = p{t) = 0. 

By considering the combination X = a + % b the main equation to solve is : 

X = -X(X 2 -Aik)/(2k 2 )] (42) 
after which we may finally obtain c thanks to ( 140]) . Solving both equations yields: 

X = (a + tb) = 2X e lt/K (iK/(X 2 e 2it/K -X 2 e~ 2it/K )) 1/2 , (43) 
cle- 2Zt/K /(X(t)X(t)) . (44) 



c 2 



Moreover the remaining transformation with A" 1 = diag(p, p, 1) allows us to fix 
k — 1 . After a final translation of the variable t, and if necessary a flip of sign we 
obtain 



a = 1/6 = y/cot(t) , (45) 



c . . (46) 

v/sin(2t) 
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Similar expressions can be obtained for Bianchi type VIIo, but they involve real 
exponentials of time. 



Bianchi class B In this case , condition (1271) is equivalent to 

1^ = 6^1,^ (e 123 = l). (47) 



For algebras of Bianchi types IV, VI^ (with h ^ 0, 1) and VII^ (with h ^ 0) we 
found that the only solution to the equations ff4T|) is I a = X a a , but for type III = 
VIi, we obtain a two-parameter solution. This is a reflection of the fact that the 
derived algebra now is of rank one instead of two. 

Equations (1241) read 

Xi pj - Xj pi = -2 a k e kij with a k = a fl 9>* k . (48) 

As we require a non singular co-frame, they are only compatible with an algebra 
of Bianchi type III. By derivation, using f|25|) we obtain the system of equations: 

Pi = {^n^g^ + 6X k p k )p t + 2 6e l jk a j p kl 

at = ^n IM/ g^ai + eX k a k p i . (49) 

To check the consistency of (1481) with (|2"U|) . let us rewrite them as relations between 
vectors of a three dimensional Euclidean space: 

2 a = -Xxp, (50) 
p = -p + 9(X- p) p + 26a x p, (51) 

Tl -* 

a = -a + 6(X-a)p. (52) 
Taking the derivative of flBTJ]) and inserting in it f lBTj [5^]) . we obtain 
2<3 + Axp = n a + 26 (X ■ a) p + — A x p + 6 (X ■ p) X x p + 26 X x (a x p), 



a + 6(X ■ p)j (2a + A x p). 



Thus if2a + Axp = 0ata period of time, it always vanishes. This insures that 
the solution obeying this condition constitutes a self-dual solution. 
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It was shown in (T3J H3] that there doesn't exist any diagonal self-dual metrics 
with a Bianchi type III symmetry group. We provide in the following the most 
general solution of type III, without any a priori (non geometrical) assumption on 
the co-frame, from which it is easy to see that diagonal metrics don't exist but 
non diagonal ones do. As previously, assuming a non singular frame, by using 
rotations and linear transformations we may put A = (1,0,0) and I a = (1,1,0). 
Then f[24"l) implies that the frame must be of the form 

(a 2c-b -2v- s \ 
m =[ r b s . (53) 

\ v c J 

This frame can still be transformed, without altering the canonical structure con- 
stants and constants of motion, by acting on the right with the GL(3, M.) matrix 

/ 1 -n n v \ 
A 1 = /i 1 - fj, -v \ (54) 
V 1 J 

and on the left with a rotation that leaves A invariant 

/l \ 

O=0 cos(0) sin(0) . (55) 
y — sin(0) cos(0) J 

The self-duality equations imply that v(t) = Vc(t). Thus we may assume that 
v(t) = after, if necessary, a rotation of angle <fi such that tan(0) = V. Moreover 
we then see that s(t) = S c(t) and that by choosing an appropriate value of v we 
may also assume s(t) = 0. The remaining parameter /i can by fixed by requiring 
6(0) = 0. 

So without loss of generality, the frame we have to consider reduces to: 

(a 2c-b \ 
9k = \ r b . (56) 




In order to pursue the integration it will be useful to redefine the unknown func- 
tions as follows (using the fact that c cannot vanish as we only consider non singular 
frames) : 

a(t) = a{t) c(t), bit) = P(t)c(t), r{t) = p(t) c{t). (57) 
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Then the self-duality equations lead to 

2 + a 2 + ap 

a = (2p-(a + P )P)c 2 ' (58) 

k = ( 5 Q) 

P (2p-(« + p)/3) C 2 ' [ * J) 

2 + p 2 + ap 

P = 7^ 7 — I — wrTT , (60) 



2 c (2 p - (a + p) f3) 2 



(61) 



They imply 



d{p - a) -( P + a) and = 4 +fr + ") 2 , (62) 



d/3 p — a 

These equations are easy to solve and lead to 

a = sinh I — =J— l — K cosh ^ 



K J \ K 

sinh I — — — | + K cosh 



Then we may express, by a quadrature, the function c as a function of (3 using the 
equation obtained from the ratio of $ and c: 

2 _ 2 COsh (ffi) 

C " C ° (/3 - 1) sinh (S) _ x cosh (S) (63) 

and finally express everything as functions of t by integrating equation (I59|) which 
leads to the surprisingly simple result (after a flip of the sign of t): 

0=*t. (64) 
Co 

The solution we discussed is generic. It depends on two arbitrary constants, namely 
Co, K and note that we have used a time translation to set /3 = at time t — 0. 
Note that we assumed that a is always different from r. If for a value of t we set 
a = r, then they become identical for all times. In this special case, the derivative 
of /3 vanishes and we obtain that b(t) = B c(t) implying the remaining equations: 

c R 
c er 

- = -A-~r («) 
a c l a 1 
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with R = ^zsy. 



To solve the equations ( |65|) we use an auxiliary function y(t) = log \a(i) c(t)\ 
and obtain, after elementary operations, the first integral: 

y 2 + AR 2 e~ 2y = L 2 , (66) 

from which we deduce the expression of the product a(t) c(t) and 

c(t) 2 = ^coth[Lt], (67) 

a(tf = -sinh[2Lt]. (68) 
lj 

Let us notice that Rt has to be positive for the solution to be real. When t = 
we encounter a curvature singularity. 



5 Self-dual Connections 

In the conventions we have used, we naturally find solutions whose curvature is 
self-dual while the connection is not. However, as we will show below, we may 
always make a local SO (A) rotation to make the connection self-dual. In fact, the 
vanishing of the anti self-dual part of the connection can be related to the inte- 
grability condition, thereby ensuring the existence of a pure self dual-connection. 
We shall illustrate this point by evaluating, in general, the gauge transformations 
which map our canonical co-frame to one whose associated connection is self-dual. 

As is well known, the two (families of) parallelisms on the three-sphere S 3 
provide a way to decompose any 50(4) rotation into the product of self-dual and 
anti self-dual SO (3) rotations. If u, x, y and z are the coordinates of a point of 
S 3 : u 2 + x 2 + y 2 + z 2 = 1, an anti self-dual SO (3) element can be written as: 
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-y 


—z 


\ 
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-y 
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—z 


u 


x 
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—x 


u 


I 



(69) 



The gauge transformation generated by such a rotation leaves invariant the self- 
dual part of the connection but transforms the anti self-dual one. Accordingly, 
when the anti self-dual part of the connection, given by the constants Ii a , is non- 
zero, such a gauge transformation may set it to zero thanks to a solution of the 
equation: 

dO = u. (70) 
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More explicitly we obtain (with := Ika& a and x l := (x,y,z)) 

du = l -x k l k , dx i = - (-u li + e ijk xj I k ) . (71) 

As a first immediate consequence we see that this gauge transformation has to be 
t independent. Moreover we also obtain from Eq. (1221) : 

dl t = l -h a c a M g & A a 1 = —4% A Ik, (72) 

which indeed ensures the integrability of the equations (171 p . In the case of a 
rank one mapping, I ka = J Q (normalized such that Ylik^i = 1); this equation 
implies there exists a function x such that d\ = I a cr a . Correspondingly, the 
gauge transformation (j69[) that leads to a frame which has vanishing anti self-dual 
connection is given by 

u = sm( X /2), x k = \ k cos(x/2). (73) 

In the case of rank three, the only possibility is Bianchi type IX. As discussed 
above it can be taken diagonal, with the matrix (Jj Q ) = — 5 ia . Moreover, it is well 
known that the components of the invariant one-forms also define parallel vector 
fields. Using the identities udu + x dx + y dy + z dz = and u 2 + x 2 + y 2 + z 2 = 1, 
it is easy to check that the invariant one-forms may be expressed as : 

a 1 = —2(xdu — udx + zdy — ydz), (74) 
a 2 = —2(ydu — zdx — udy + xdz), (75) 
a 3 = —2(zdu + ydx — xdy — udz), (76) 

which satisfy the invariance relations: da 1 = £ % jk a ^ A c fc as well as (j71|) . In the case 
where the matrix (J iQ ,) is not diagonal, it must be (in order to satisfy to equations 
([2"6"1) ) of the form (7^) = — (Oj Q ) with O an 50(3) rotation matrix. It is elementary 
to verify that the equations (I7ip are solved by similar functions involving variables 
u and (x fk ) = (x',y',z') related to (x k ) = (x,y,z) by the rotation O: 

x'k = -I M x l (77) 

as it must be. 

6 Conclusion 

In this note we have achieved a complete classification of all the self-dual Euclidean 
spaces admitting a G3 simply transitive homogeneity group. For each Bianchi 
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group, there are as many classes of solutions as homomorphisms I : Q3 — > so (3). 
For each of these homomorphisms, labeled by its rank, the self-duality constraint 
is described by a specific set of equations - like Lagrange or Darboux-Halphen 
in Bianchi IX, leading to distinct solutions such as Eguchi-Hanson or Taub-NUT 
respectively. These equations have non-everywhere vanishing solutions for all class 
A Bianchi groups (I, II, VIo, VIIo, VIII and IX), whereas in class B solutions ex- 
ist only for Bianchi III. The latter solutions (Eqs. flHT]) . f )83|) ) are all rank-1 and 
non-diagonal, contrary to the lore that it seamed "very unlikely to construct non- 
diagonal self-dual solutions of Bianchi types I-IX" [13] H%] . This expectation was 
empty for Bianchi A because in the present formulation, the diagonal ansatz is the 
most general, but not for Bianchi B. 

It is fair to stress that the above conclusions have been reached by adopting a 
metric representation, Eq. (JTJ), such that the foliation of Ai as K. x S3 is manifest 
and adapted to the splitting of the group SO (A) into self-dual and anti-self-dual 
factors. Other choices for the metric may exist, where the distinction between the 
various classes of self-dual solutions is less sharp. This happens e.g. for Bianchi 
IX, in the Gibbons-Hawking representation [19], where Eguchi-Hanson (rank-0) 
and Taub-NUT (rank-1) are indeed, to some extent, unified. This is possible, 
however, at the price of abandoning the G3-invariant frames that we use. Here, 
this choice has been instrumental, not only for providing the classification in terms 
of the rank of the g3 — » so (3) homomorphism, but also for scanning all possible 
solutions. 

Our analysis can easily be adapted to the ultra-hyperbolic case (spacetime sig- 
nature (— , — , +, +)), where the algebra of Bianchi type VIII will play a role anal- 
ogous to the one of Bianchi IX here. It should however be emphasized that we will 
have to distinguish, in that case, among the inequivalent one- dimensional subalge- 
bras of so(2, 1). Other extensions concern the addition of a cosmological constant 
(Weyl self-duality), or the implementation of the method to higher- dimensional 
set-ups admitting self-duality, as e.g. in seven dimensions. Last but not least, the 
physical analysis of the new Bianchi III solutions remains to be completed. The 
large number of moduli in flHT]), (183]) make such an analysis quite involved. 

It would be interesting to discover if the techniques and classification presented 
in this paper can be generalized to cases where we have a more general setup than 
vacuum Einstein equations. For example, the evolution of certain Bianchi spaces in 
the presence of a cosmological constant, as well as various matter such as a perfect 
fluid and/or scalar, spinor, and electro-magnetic fields has been studied recently 
in [22] [23] E3J. The authors have derived exact solutions in certain cases, some 
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of which are singularity free, and a complete classification of solutions would be 
appreciated. In the case of Einstein's equations with a cosmological constant, it is 
well known that requiring self-duality of the Weyl tensor is sufficient to satisfy the 
equations of motion. However in this case our technics have to be abandoned as 
we are generically confronted with a system of second order differential equations. 

We note that most of the Bianchi classes have solutions which necessarily have 
singularities. It would be interesting to discover if these singularities could be cured 
by the introduction of either a cosmological constant and/or additional matter 
fields. Examples have already been discussed in [221 [231 121] as well as [25l [26] . 
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Appendix: G3 self-dual metrics 

We summarize here all known (real) self-dual G3-homogeneous gravitational in- 
stantons, following the classification pattern we have developed in terms of ho- 
momorphisms 03 — > so (3). These can be of rank 0,1 or 3. Whenever rank-0 
and rank-1 solutions coexist (i.e. for Bianchi I, II, VIo and VIIo), they are both 
captured by a single expression with a two- valued parameter e = 0,1. Bianchi 
IX is the only case which possess rank-0 and rank-3 solutions. Bianchi III is the 
most peculiar: contrary to the other classes, it requires a non-diagonal co-frame 
(T3TT) (and consequently a non-diagonal metric ([1])) and allows exclusively rank-1 
solutions. These solutions are new. 

In the following presentation, integration constants that can be reabsorbed by 
coordinate redefinitions have been discarded. Consequently, all remaining param- 
eters are genuine moduli of the solutions. 
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Bianchi I The solution reads 



ds l = e" 2£i dt 2 + e" 2et (a 1 ) 2 + (a 2 ) 2 + (a s ) 2 (78) 

and describes flat geometry for both rank-e solutions (e = 0, 1). Notice, however, 
that the rank-0 case appears naturally in Cartesian coordinates, whereas rank-1 
emerges in a kind of mixed Cartesian/polar, Euclidean-Rindler-like coordinates. 

Bianchi II The rank-e solutions are given by 



ds 2 = '-^dt 2 + V) 2 + Ua 2 ) 2 + ^-^(a 3 ) 2 . (79) 



The Kretschmann scalar is given by 



K = Rmnp Q R MNPQ = 8b 2 4e- 4t ^ {et)2 + 3b t ° et + 3bl 
from which we can see that at t — the metric has a curvature singularity. 



50) 



Bianchi III The Bianchi-III self-dual metrics are all rank-one. They are cap- 
tured by two expressions, depending on whether a ^ r (Eq. flHTl) ). or a = r (Eq. 
f[83l ). in the general co- frame fl31|) after a convenient rescaling of t by t — )■ c^t we 
obtain: 

dg2 = F(t) df2 (Kt-l)tanh(t-U)-K 

32 c 4 cosh 4 (t-t t ) c 2 

,9l{t) f 1X2 , 92{t) ( 2N2 , „ f3W I 2 /o, \ 

where we have introduced the functions -F(t), <7i(£)> 92(f) and (73 (i) given by 

F(t) = 8 c 2 cosh (t - t*) (A' i - l) sinh (t - t*) - K cosh (t - t*) 

^(t) = (AT 2 + l)cosh(2(t-t»)) +2ATsinh(2(t-t*)) +2A 2 t 2 + A 2 - 1, 
£ 2 (t) = (A 2 + 1) cosh (2(t - Q) - 2K sinh (2(t - Q) 

+2K 2 1 2 -8Kt + K 2 + 7, 
g 3 (t) = (K 2 -1) cosh (2(f - +2 A 2 t 2 -AKt + K 2 + 1. (82) 

Denoting the zeroes of F(£) as ti, we find that the Kretschmann scalar diverges as 
(t — ti) 6 as t — > U, indicating a curvature singularity at each of the two zeros of 
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F(t). The second solution with a = r is given by 

ds 2 = L 1} tanh 2 ( ^ N dt 2 + L(B-1) tanh(t) (a 3 ) 2 
4 cosh 2 (t) v ' v 1 

+ 4(/-l) {[ 2i?2csch (2t) + tanh(t)] (a 1 )^ 

-2 [tanh(t) — 2 B (B — 2)csch(2 *)] cr 1 a 2 

[csch(2 t) (cosh(2 1) + 2 (5 - 2) 2 - l] (a 2 ) 2 } . (83) 

In this case, the Kretschmann scalar is given by 

_ 384 coth 6 (t) 

K ~ (B-iyL* (84) 

and we find at t = the metric has a curvature singularity. The general analysis 
of these geometries is involved and provides interesting features, which deserve a 
separate study 

Bianchi VIo The general metric for rank-e solutions reads: 

ds 2 = sin(2t) < —^- (dt 2 + (a 3 ) 2 ) + tan(t) (a 1 ) 2 + cot(t) (a 2 ) 2 , (85) 

where t G [0,7r/2]. A curvature singularity appears at the boundaries t = and 
t = 7r/2. The Kretschmann scalar is given by 



, 4 _ 4feZ 6eZ sin(4t) + [9 + (eZ) 2 ] cos(4t) + (e Z) 2 + 15 

/C = 16c e ■ 6/ ^ ' ^ 

sin (2 t) 

and we see that at t = 0, tt /2 the metric has a curvature singularity. 

Bianchi VIIo The general metric for rank-e solutions reads: 

ds 2 = sinh(2t) (dt 2 + (a 3 ) 2 ) + tanh(t) (a 1 ) 2 + coth(t) (a 2 ) 2 (87) 



c6 



The Kretschmann scalar is given by 



/C = 16 c n e 



UeZ 6eZ sinh(4t) + [9 - (e Z) 2 } cosh(4t) - {eZf + 15 

^ ■ 6~~ 



sinh b (2£) v ' 

and we again find a curvature singularity at t — 0. In the rank-1 case, we may 
take the large-t limit of the above metric to obtain 

ds 2 = cle 2 ^ 1 [dt 2 + (a 3 ) 2 ) + (a 1 ) 2 + (a 2 ) 2 (89) 

which also solves the self-duality equations but has a vanishing Kretschmann scalar 
and is actually flat. 
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Bianchi VIII There is only the rank-0 solution whose simplest writing is: 

, 2 P~ 1 ' 2 ,2 1/9 ( (^) 2 (^) 2 \ , ^ 

ds 2 = ——dx 2 + P 1 ' 2 + + 90 

4 \xi — x x 2 — x x J 

with P = (xi — x)(x 2 — x)x where x±, x 2 are positive and < x < min{xi, x 2 }. 
The Krestschmann scalar blows up as x~ 3 (x — xi)~ 3 (a; — x 2 )~ 3 , indicating that 
the singularity at the boundary of the domain of definition of x is a true one. 
Let us assume X\ = a 2 > x 2 = j3 2 . We may re-express the metric in the gauge 
N = 6, using Jacobi elliptic functions of module (3 /a. Making the substitution 
x = (3 2 sn 2 (at) we obtain: 

ds 2 = a(3 2 sn(at) cn(at) dn(at)dt 2 (91) 

(3 2 sn(ai) cn(at) 1 2 asn(at) dn(at) 2 2 acn(at) dn(at) 3 2 

H ~j 7 7\ °" ) H 1 T\ \ a ) H 7 it (o - ) • 

adn(at) cn(at) sn(at) 

The t variable belongs to the interval [0, 2 K(f3 / a) / a], K{k) being the complete 
elliptic integral of the first kind of module k. 

Bianchi IX We have solutions in case of rank-0 (Lagrange system of equations) 
and of rank-3 (Darboux-Halphen system of equations). 

The general rank-0 solution was found by Belinsky et al in [5] as a strict 
S'[/(2)-symmetric generalization of the Eguchi-Hanson gravitational instanton. 
The latter has enhanced SU(2) x U(l) isometry and is a solution of the Lagrange 
system. An algebraic expression of the solution, analogous in form to Eq. (1901) . is 

4 \x — xi x — x 2 x J 

with P = (x — Xi)(x — x 2 )x where Xi, x 2 are negative. The general strict SU(2)- 
symmetric solutions have curvature singularities. To obtain the Eguchi-Hanson 
gravitational instanton, one may take x\ = x 2 = xo and by a translation in x set 

X 3 = 0. 

Let us assume x% = —a 2 < x 2 = —f3 2 . In the gauge N = 9, the metric can 
still be expressed in terms of standard Jacobi elliptic functions but now of module 
\/a 2 — f3 2 /a; using x = (3 2 sn 2 (a t) /cn 2 (at) with t G [0, K(^Ja 2 — f3 2 /a)/a] we 
obtain: 

ds 2 = a ^ n{at ]^[ at) dt 2 (93) 
cn 3 (ai) v ; 

(3 2 sn(at) 12 a sn(a t) dn(a t) 22 adn(at) 32 

H 7 r— : — 7 r((T ) -\ (a ) H r r-((T ) . 

acn(ai) dn(at) cn(at) sn(at)cn(at) 



7 The metric can be chosen diagonal; the off-diagonal entries are consistently set to zero in 
this class. 
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The issue of the rank-3 solution is more subtle. It is a generalization of Taub- 
NUT metric, with strict SU(2) isometry, and solves the Darboux-Halphen system. 
As explained in [8], it requires the use of modular forms. For simplicity we trade 
a,b,c for Q a , a = 1,2,3, where Qi = Q 2 = — , and f2 3 = -k-, and introduce a 
triplet of weight-two modular forms of T(2) C PSL(2, Z): 

dX/dz dX/dz dX/dz 

El = T ' E2 = A^T Es = A(A = 1) ' (94) 

where A is solution of Schwartz' equation 

d 3 X/dz 3 3 f d 2 X/dz 2 \ 2 _ 1 1 \ UXV 

dX/dz ~ 2 V dX/dz ) ~ ~ 2 V A 1 + (A - l) 2 ~ A(A -I)) \~dz~ ) ' ^ ' 

Any real solution of the Darboux-Halphen system reads: 

n a (t) = -~]ogE a (it), a = 1,2,3, (96) 

and generates a Bianchi-IX, rank-1 gravitational instanton (Eqs. ([T]) and (1311) ). A 
generic curvature singularity is present in these geometries. This singularity can 
be pushed at infinity for one specific choice of X(z): X = Q\jQ\^ where #2,3 are the 
standard Jacobi theta functions^). This case corresponds to the Atiyah-Hitchin 
solution (201 EE] • Similar to the rank-0 solution, in the case that two of the metric 
factors are equal the solutions has enhanced 577(2) xU(l) isometry and takes the 
form of the well-known Taub-NUT metric 

" = + <* 2 > 2 ) + j^ 2 + irrW* 2 (97) 
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